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Abstract 

Dimensional reduction occurs when the critical behavior of one system can be re- 
lated to that of another system in a lower dimension. We show that this occurs for di- 
rected branched polymers (DBP) by giving an exact relationship between DBP models in 
D + 1 dimensions and repulsive gases at negative activity in D dimensions. This implies 
relations between exponents of the two models: ^{D + 1) = a{D) (the exponent de- 
scribing the singularity of the pressure), and u±{D + 1) = i^{D) (the correlation length 
exponent of the repulsive gas). It also leads to the relation 9{D + 1) = 1 + (t{D), where 
(t{D) is the Yang-Lee edge exponent. We derive exact expressions for the number of 
DBP of size N in two dimensions. 

PACS numbers: 64.60.Fr, 04.20.Jb, 04.60.Nc, 05.20.Jj 

The phenomenon of dimensional reduction has attracted considerable attention over the years. 
The first example was the controversial random field Ising model (RFIM), whose critical 
behavior was conjectured to be the same as the pure Ising model in two fewer dimensions | ll. 
A proof of long-range order for the RFIM in three dimensions ||2l|3]| showed that dimensional 
reduction fails there, and recent work llH |5l |6l| has elucidated what goes wrong. A second 
example is the Parisi-Sourlas reduction of branched polymers (BP) m D + 2 dimensions 
to the Yang-Lee edge or iip'^ field theory in D dimensions [7|. This was recently confirmed 
with the discovery of an exact relationship between BP models and repulsive gases at negative 
activity in two fewer dimensions |l8l|9||. The failure of the heuristic arguments for dimensional 
reduction in the RFIM underscores the importance of having an exact result. In this letter we 
give a third example, in which directed branched polymers (DBP) reduce to the repulsive gas 
at negative activity in one fewer dimension. 

We consider directed branched polymers as self-avoiding tree graphs embedded in Z^+^ 
or MP^^ so that every vertex can be reached from the root at by a sequence of links which 
move forward with respect to a preferred direction. See Figure [T] Let denote the number 
of DBP with N vertices, and let Zdbp(-2) = Y1,n dNZ^ ■ We prove the identity 

Pnc{z) = -Zdbp(-^), (1) 
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Figure 1 : A directed branched polymer in 1? . 



which relates the DBP generating function to the density of an associated repulsive gas model. 
In contrast to the undirected case, the reduction in dimension is one. If we define a from the 
singularity pHC ^ [z — -2c)^~" of the hard-core gas at the negative activity critical point, and 
7 from the singularity Zdbp ~ (2; — 5c)^^'^, then we find = —Zc and 'j{D + 1) = a{D). 
Define an exponent 9 from the asymptotic behavior ~ z~^N~^. Then 6* = 2 — 7. Noting 
that the Yang-Lee edge exponent a Ifim can be identified with 1 — a lITTlfT^ . we obtain the 
relation 

e{D + l) = l+a{D). (2) 

In one and two dimensions, a can be determined from the exact values a{l) = | ifTTl . 
a (2) = I (the latter follows from the solution to the hard-hexagon model [13| at the nega- 
tive activity critical point, see [.14t.l5J : alternatively from [16], assuming the model is in the 
Yang-Lee class.) Hence ^^(2) = | and ^(3) = |. We also obtain identities relating DBP corre- 
lations with repulsive gas correlations, which imply that the DBP exponent for the transverse 
correlation length ^±{0 + 1) equals the repulsive gas exponent i^{D). 

The related problem of directed lattice animals (DA) has been studied extensively ifTTl 
[ia[I3|20linil23[Il[I3|23l2l|23|26l|27J. It is generally believed that the loop-free 
condition does not affect the critical behavior, so that both systems should have the same 
exponents. For DA, there are exact results in two dimensions (see ll22l l23l l24l . the review 
ll25i . and references therein) and in three dimensions [14|. Also, in any dimension, models 
of DA have been related to dynamical models of hard-core lattice gases (14] (see also Ii28i . 
which connects Lorentzian semi-random lattices with DA and gives an alternative derivation 
of Dhar's equivalence). They have also been related to the critical dynamics of the Ising 
model in an imaginary field ifTTlfTSl . As a result, the identity Q is believed to hold for DA 
EIEIIISI, and the values in D = 2, 3 should be the same as the ones given above for DBP. 

We consider a class of DBP models where the vertices yk = {tk, Xk) have a time compo- 
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nent E M+ and a space component G S, where S is either or Z^. If a lattice model 
is desired, the time component can be discretized by taking limits (see examples below). Let 
T be a tree graph on {1, ... , A^}, and let yk be the position of the k^^ vertex. Fix the vertex 1 
as the root, with yi = 0. For each pair (i, j) = ij, define 

Each link of T connects a vertex j to a vertex i, where i is one step closer than j to the root 
along T. As we are considering directed BP, we require tj > ti. 

The weight associated with each DBP configuration depends on a linking weight V{y) = 
V{t,x) and a repulsive weight U{y) = U{t,x) = U{t, —x). The generating function for 
(rooted) DBP is written as 

oo AT p 

Z^bA^) = J2 7]v^ E / n ldyJ^V{y,,)] n U{y,.). (4) 

Here each pair {i, j} appears exactly once, either in HjieT °^ Hji^T- assume U{t, x) 
1 if either t or x tends to infinity (repulsion vanishes at infinite distance). In order to get 
dimensional reduction, we require 

V{t,x) = U\t,x), (5) 

where prime denotes the t-derivative. For positive weights, we require that U and U' = V 
are positive. Also, V needs to be an integrable function of y for Q to be well-defined. 

Our main result relates Zdbp to the density of a repulsive gas in the /^-dimensional space 
S. Let A C S'. Using the grand canonical partition function 

oo „ N 

Z^c{z) = Y,'— \{dx, n U{Q,x,,), (6) 

Ar=0 ■ "^^^ i=l l<j<7<Ar 



we can define the density 



Pnc{z) = z^- 



lim -^logZnciz) 



(7) 



dz [a/s |A| 

Theorem 1. For all z such that the right-hand side converges absolutely, 

Prc{z) = -Zbbp{-z). (8) 

We now give some continuous and discrete examples of models of type ©. If one takes 
U{t, x) = ^(t + |a;| — 1), where i!) is the usual step function, then V{t, x) = 6{t + ^\ — 1). 
One obtains a DBP model where each monomer is a hard diamond (when |a;| = X]q=i l^al) 
or a hard double-cone (when |xp = J2a=i^Q)^ obtained as {{t,x) : \t\ + \x\ < 1}. Each 
monomer is distributed uniformly in contact with the positive surface of the monomer it 
is linked to (subject to the constraint of nonoverlap with other monomers). One can also 
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consider a hard-sphere model by taking U (t, x) = diiP' + |a;p — 1). Then the monomer is still 
distributed uniformly in x, which means angles close to the preferred direction are favored. 

Theorem [T] equates the generating functions of these {D + l)-dimensional DBP models 
to the density of the hard-sphere gas in D dimensions, at negative activity. For D = 1, the 
pressure of the hard-rod gas is computable; it is 

p{z) = LainbertW(2) = -T{-z), (9) 

where T{z) = Y1'n=i N^~^ /N\ is the tree generating function Q. Hence 

Zbbp{z) = -PHC(-^) = Yl ^TT- ^^^^ 

N=l 

This gives a simple expression for the volume available to DBP of size A^; this can be checked 
with some effort for small values of A^. 

A natural class of lattice examples can be obtained by taking U{t,x) = 1 — I{x)'d{l — t), 
where I{x) is an indicator function of a set of "neighbors" in the lattice, such as {x : |x| < 1}. 
Since V{t,x) = I{x)5{t — 1), the set determines which sites a link can jump to; t always 
increases by 1. Figure [T] shows a representative DBP for this model in two dimensions. On 
the gas side of identity ([T]), t is set to 0, and we obtain various nearest-neighbor exclusion 
models, including the hard-hexagon model in D = 2. When counting allowed DBP, it is 
important to include the factors Hji^T ^ which enforce the nearest-neighbor exclusion 
for monomers on the same level (same value of t). Also, if a monomer at level t has n 
neighbors at level t — 1 in the polymer, the associated [/-factors can be written as t9"~^(0), 
which should be interpreted as ^ = J (as can be seen by approximating with a 

smooth function and integrating over t). Monomers in levels more distant than this do not 
interact (similar conditions occur in Theorem III.2 of 129.1 : in fact (fT3t can be used to give an 
alternate proof of that identity). 

Theorem [U equates the generating functions of these models with the nearest-neighbor 
exclusion models in D dimensions associated with the weights f/(0, x) = 1 — I{x). For the 
hard-hexagon model in D = 2, there is an exact solution for the pressure llT3l . In D = 1 we 
have a dimer model, for which the pressure is known, see [11, Equation 2.16]: 

p(2) = In (I + |yrT4i) . (11) 

Thus the generating function is 

^ . ^ d , , ^ [2iV-l]!!2^-i^^ 

N=l 

which enumerates the number of DBP with monomers. 

In order to prove Theorem [T] we derive an identity which encapsulates the effect of inter- 
polating in the t variables. Let /(t) be a smooth function of {tj}, {tjj} which approaches 
when any of the t's tends to infinity. Then di. forest-root formula holds: 

/(o) = Y. L n [-^('^^^ - ^^)]/^^'''Ht). (13) 

{F,R) '^*-+ reR jieF 
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Here R is called the set of roots and is any subset of {1, ... , A^}. We sum F over the set of all 
forests or loop-free graphs on {1, . . . , A^} with the property that each connected component 
or tree of F contains exactly one root. Each link in a tree of F connects a vertex j to a vertex 
i which is one step closer to the root for that tree. The integration region is {tr > 0, r E 
R and tj > ti, ji E F}. 

For = 1, (US) reduces to /(O) = - f'{t)dt. For N = 2, consider /(ti, ta, ^12) and 
use subscripts 1, 2, 12 to denote partial derivatives. Then 

POO 

/(O) = - / dsihis, s, 0) + /2(s, s, 0)). (14) 

Apply the = 1 formula to the /i term, integrating with respect to X2 — xi, and to the /a 
term, integrating with respect to xi — X2. The result is 

/•oo /"OO 

/(O) = / dtj d{t2 - tl)(/l,2 + /l,12) + (1 — - 2), (15) 

Jo Jo 

since = 1 for t2 > h, and = 1 for ti > ^2- The two /i_2 terms combine to form 
/jg2 dtidt2fi^2, which is the term R = {1, 2} of ([T3t . The other two integrals are the terms 
R = {1},{2}. 

One can prove the general case by induction on A^. Begin as above with 

/(0) = - / dsJ2Ms,...,s,0,...,0), (16) 
-^0 k=l 

where the integral is along the diagonal, ti = t2 = ■ ■ ■ = tj^. Consider one of these terms, 
say k = N, and apply (fT3t in the variables = — ^at, z = 1, . . . , A^ — 1, keeping = s 
fixed: 

MtN,...,t^,0,...,0)= ^ f^^U^-dir] ll[~d{i,-ii)]fif^^\t). (17) 

{F,R) ^+ rdk jidF 

Note that when computing the derivative of with respect to there will be a term /at and 
also a term fN,rN (coming from the dependence on trN = U — tjy = ir). Thus each (F, R) 
on {1, . . . , A^ — 1} gives rise to 21'^' terms, each of which can be assigned a unique (F, R) 
on {1, ... , A^}. R consists of A^, together with each r E R with an fjy^r term. F consists of 
F, together with rN, r E R when r gives rise to an fN,rN term. Observe that each root in R 
ceases to be a root in R if it is connected by a bond rN in F. We obtain in this way all (F, R) 
with N E R, and each satisfies the condition that each tree of F contains exactly one root. 
As a result, 

N 

/(0) = E E J-dt,] n [-d{U-h)]l[[-d{t,-U)]f(^'^\t). (18) 

k=l {F,R):keR ^+ i€R\{k} ji&F 

It is evident that if we take the term (F, R) of ( fT?l . and consider the subset of the integration 
region for which = min^-g^ tr, we obtain the term k, (F, R) of (fTSl) . 
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The forest-root formula (fOl) is the key to dimensional reduction; in fact a two-dimensional 
version of (fT^ was used in [8, 9| to derive a repulsive gas mapping for ordinary (isotropic) 
BP. The argument proceeds by applying ([T3t to 

N 

fit)=g{h/e)l[gieU) J] U{U„Xi,), (19) 

1=2 l<i<j<N 

where g is any smooth function which decreases to and satisfies (7(0) = 1. Thus, 
PHcW = limlim^-^X;- j n^x./(0) 

= ^ z^)^wX ni-''!'.! n ^<^>^n- (20) 

For each G F, U{yji) is differentiated and becomes the linking weight V{yji). For each 
r G -R, a is differentiated: 

/(^'«)(t) = n^'H^ n ^(^:'-^) n ^(%-^)- 

Note that the distribution of for an n-vertex tree is essentially —{g{etrY)'dtr, which is a 
very spread-out probability measure. One finds that due to the large separation in the t di- 
rection, all the interaction terms U between different trees are equal to 1 (up to terms which 
vanish as e ^ 0). This means that the sums/integrals associated with each tree are inde- 
pendent of the others, and all the trees of F decouple. One tree has its root fixed at by 
a factor —{g(ti/e))', which converges to 5{ti). The others cancel with the normalization 
Znc{z), so that (l20b reduces to Q, with an additional —1 for each link of the tree. Thus 
Phc(-2) = —Zdbp{—z), which is Theorem[ll 

One can prove identities relating n-point functions of the repulsive gas to n-point func- 
tions of DBP. If one takes the derivative of (U]) with respect to a source at x, one obtains for 
the 2-point correlation functions 

Ghc(0,x; z) = - / Gdbp(0,?/; -^)(it, (22) 
Jo 

where y = {t, x) G M+ x M^. Hence the transverse correlation length exponent z/x(-D + 1) 
of the repulsive gas must be the same as i^{D). The value z/(l) = ^ can be computed directly, 
and the value z/(2) = follows from hyperscaling (Dv = 2 — a) (or equivalently from 
Dh'± = 2 — a QUI Equation 28]) with a = | as indicated above. One can also consider 
unrooted DBP (divide each term in @ by N) and relate them to the pressure of the associated 
repulsive gas. 

In conclusion, we have demonstrated the underlying mechanism of dimensional reduction 
for directed branched polymers. In the process, some lattice and continuum models of DBP 
are solved exactly in two and three dimensions by reference to repulsive gases in one lower 
dimension. 
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